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ABSTRACT 

We derive a new and user friendly simple analytic approximation for deter- 
mining the effect of intergalactic absorption on the 7-ray spectra of TeV blazars 
the energy range 0.2 TeV < E 1 < 2 TeV and the redshift range 0.05 < z < 0.4. 
In these ranges, the form of the absorption coefficient t(E 7 ) is approximately 
logarithmic. The effect of this energy dependence is to steepen intrinsic source 
spectra such that a source with an approximate power-law spectral index T s is 
converted to one with an observed spectral index r o ~ T s + AT(z) where AT(z) 
is a linear function of z in the redshift range 0.05 - 0.4. We apply this approxi- 
mation to the spectra of 7 TeV blazars. 

Subject headings: Gamma-rays: general - blazars: individual (Mkn 180, 
PKS 2005-489, PKS 2155-304, H 2356-309, 1ES 1218+30, 1ES 1101-232, PG 
1553+113) 

1. Introduction 

Stecker, Malkan & Scully (2006) (SMS) have used recent Spitzer observations (Le 
Floc'h et al. 2005, Perez- Gonzalez et al. 2005) along with other data on galaxy luminosity 
functions and redshift evolution in order to make a detailed evaluation of the intergalactic 
photon density as a function of both energy and redshift for < z < 6 for photon energies 
from .003 eV to the Lyman limit cutoff at 13.6 eV in a ACDM universe with Q A = 0.7 
and il m = 0.3. They then used their calculated intergalactic photon densities to calculate 
the optical depth of the universe, r, for 7-rays having energies from 4 GeV to 100 TeV 
emitted by sources at redshifts from to 5. They also gave a parametric fit with numerical 
coefficients for t{E 1 , z). 
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As an example of the application of the detailed numerical determination of the optical 
depth, they calculated the absorbed spectrum of the blazar PKS 2155-304 at z — 0.117 
and compared it with the spectrum observed by the H.E.S.S. air Cherenkov 7-ray telescope 
array. It was noted that a steepening in an E~ 2 power-law differential photon source 
spectrum for PKS 2155-304 would be steepened by approximately one power to E~ 3 as 
would be produced by a 7-ray opacity with a logarithmic energy dependence. 

The purpose of this letter is to generalize this result by determining fits for 
approximating t(E j ,z) by logarithmic functions in order to predict power-law steepenings 
in assumed blazar source spectra of the form E~ Ts to observed spectra of the form E~ r ° 
where T ~ T s + AT and Ar is determined to be a linear function of z. 



2. Calculation 

In order to derive our results, we fit the results from SMS for r(_E 7 , z) to a form which 
is assumed to be logarithmic in in the energy range 0.2 TeV < E 7 < 2 TeV and which 
has a linear dependence on z over the range 0.05 < z < 0.4. It is important to note that 
our linear fit to the z dependence is both qualitatively and quantitatively different from the 
linear dependence on redshift which would be obtained for small redshifts z « 1 and which 
simply comes from the fact that for small z where luminosity evolution is unimportant and 
where r oc d, with the distance d ~ cz/H oc z, our quantitaitve fit for the higher redshift 
range 0.05 < z < 0.4 comes from the more complex calculations based on the models of 
SMS. For this reason, our linear fits are not simply proportional to redshift. 

We choose the form of our fitting function inn both redshift and energy to be given by 



t(£ 7 , z) = (A + Bz) + (C + Dz) ln[£ 7 /(l TeV)], 
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where A, B, CandD are constants. 

If we then postulate an intrinsic source spectrum which can be approximated by a 
power law over this limited energy range of one decade, 

the spectrum which will be observed at the Earth following intergalactic absorption will be 
of the power-law form 

$ (£ 7 ) = Ke- iA+Bz ^E-^ +c+Dz \ 

This can be compared with the empirically observed spectrum which is usually 
presented in the literature to be of the power-law form. The observed spectral index, T Q , 
will then be given by 

r = r s + at(z) 

where the intrinsic spectral index of the source is steepened by AT(z) = C + Dz. 



Table 1. Optical Depth Parameters 



Evolution Model 


A 


B 


C 


D 


Fast Evolution 


-0.475 


21.6 


-0.0972 


10.6 


Baseline 


-0.346 


16.3 


-0.0675 


7.99 
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The parameters A, B, C, and D obtained by fitting the optical depths derived for the 
fast evolution (FE) and baseline (B) models of SMS are given in Table 1. Figures [T] and [2] 
show the excellent fits of these parameters to a linear dependence in both In .EL and redshift. 
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Fig. 1. — Linear fits in lnE 1 to the optical depths calculated in SMS (points shown) for 
their fast evolution (FE) case in the energy range 0.2 TeV to 2.0 TeV for redshifts of (bottom 
to top) 0.05, 0.10, 0.15, 0.20, 0.25, 0.30, 0.35, and 0.40. A similar result is obtained for the 
SMS baseline (B) case. 



3. Spectral Indices of Individual Sources 

Table 2 gives a list of blazars which have been detected at TeV energies and for which 
spectral indices (T ) have been measured in the energy range 0.2 - 2 TeV. We also give 
the observed redshifts of these sources and the intrinsic spectral indices of the sources (r s ) 
derived from the baseline (B)and fast evolution (FE) models using our analytic expressions 
for AT(z). A spectral index of less than 2 indicates that the energy range of the observation 
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Fig. 2. — The fits obtained for the linear functions A + Bz and C + Dz shown for the FE 
case as descibed in the text. A similar result is obtained for the B case. 

would be below the Compton peak energy in the spectral energy distribution, E^(E 1 ) of a 
synchrotron self-Compton (SSC) model (Stecker, De Jager & Salamon 1996). In the case of 
the blazar PG 1553+113, the energy range of the observation is from 0.09 to 0.4 TeV which 
is somewhat below the energy range we have used in deriving our approximation. 



We have derived a new and user friendly simple analytic approximation for determining 
the effect of intergalactic absorption on the spectra of 7-ray sources in the energy range 0.2 
TeV < E 7 < 2 TeV and the redshift range 0.05 < z < 0.4, where the absorption is primarily 
from interactions with intergalactic photons in the optical to near infrared wavelength 
range. In these energy and redshift ranges, the form of the absorption coefficient t(E 7 ) is 
approximately logarithmic. The effect of this energy dependence is to steepen the intrinsic 



4. 



Conclusions 
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Table 2. Spectral Indices for Seven Observed Blazars in the 0.2 - 2 TeV Energy Range 



Source 


z 


r, 


* 




r, (B) 


r s (fe) 


Mkn 180 


0.045 


3.3 


(a) 


3.0 


2.9 


PKS 2005-489 


0.071 


4.0 


(b) 


3.5 


3.4 


PKS 2155-304 


0.117 


3.3 


(c) 


2.4 


2.2 


H 2356-309 


0.165 


3.1 


(d) 


1.9 


1.5 


1ES 1218+30 


0.182 


3.0 


(e) 


1.6 


1.2 


1ES 1101-232 


0.186 


2.9 


(f) 


1.5 


1.0 


PG 1553+113 


0.36* 


4.2 


(g) 


1.4 


0.5 



*MAGIC and H.E.S.S. references for observed spec- 
tral indices denoted by letters are (a) Albert et al. 
(2006b), (b) Aharonian et al. (2005b), (c) Aharonian 
et al. (2005a), (d) Aharonian et al. (2006a), (e) Albert 
et al. (2006a), (f) Aharonian et al. (2006b), (g) Albert 
et al. (2006b). 

^Although Albert et al. (2006b) and Aharonian et al. 
(2006c) use assumptions about the source spectrum of 
this object to place upper limits on z of 0.78 and 0.74 
respectively, the catalogued redshift of PG 1553+113 is 
0.36 (Hewitt and Burbidge 1989). 
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source spectrum such that a source with an approximate power-law spectral index T s 
is converted to one with an observed spectral index r o ~ T s + Ar in the energy range 
0.2 - 2 TeV, where AT(z) is a linear function of z in the redshift range 0.05 - 0.4. We 
have applied this approximation to the spectra of seven TeV blazars. These power-law 
approximations, both observational and theoretical, are only useful over the limited energy 
range indicated, viz. 0.2-2 TeV. The actual spectra should exhibit some curvature. At 
higher energies the spectra should cut off either because of a natural upper limit to the 
source spectra or, more likely in the case of high frequency peaked BL Lac objects (HBLs), 
because of the increase in optical depth with energy owing to interactions with the much 
more numerous intergalactic far infrared photons. A full numerical treatment of absorption 
can be made using the exact form of the optical depth as a function of redshift and energy 
as calculated in the detailed work of SMS. Examples of more detailed treatments of effects 
of intergalactic 7-ray absorption are given for Mkn 501 and Mkn 421 by Konopelko et 
al. (2003) and for PKS 2155-304 by SMS. However, because of the limited energy range 
of empirical TeV spectral determinations, and because of the significant statistical and 
systematic uncertainties in the empirical determinations of spectral indices of TeV sources, 
our approximate linear relations can be quite useful in making rapid and simple analyses of 
TeV source spectra. 
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